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Abstract 

Motivated by several observational and theoretical developments concerning the 
variability of Newton's gravitational constant with time G(t), we calculate the 
varying G correction to the statefinder parameters for four models of dark energy 
namely interacting dark energy holographic dark energy, new-agegraphic dark 
energy and generalized Chaplygin gas. 
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I. INTRODUCTION 



Recent cosmological observations obtained by SNe la WMAP [2|, SDSS [3| and X-ray 
{4] indicate that the observable universe experiences an accelerated expansion. Although 
the simplest way to explain this behavior is the consideration of a cosmological constant p], 
the known fine-tuning problem 6| led to the dark energy paradigm. The dynamical nature 
of dark energy, at least in an effective level, can originate from a variable cosmological 
"constant" [7J, or from various fields, such is a canonical scalar field (quintessence) |8j, a 
phantom field, that is a scalar field with a negative sign of the kinetic term 9[, or the 
combination of quintessence and phantom in a unified model named quintom [lO]. Finally, 
an interesting attempt to probe the nature of dark energy according to some basic quantum 



gravitational principles is the holographic dark energy paradigm 

Dirac's Large Numbers Hypothesis (LNH) is the origin of many theoretical studies of 
time-™ G. According to LNH, the value of G/G should be approximately the Hubble 
rate [121 ] . Although it has become clear in recent decades that the Hubble rate is too high to 
be compatible with experiments, the enduring legacy of Dirac's bold stroke is the acceptance 
by modern theories of non-zero values of G/G as being potentially consistent with physical 
reality. Moreover, it is the task of a final quantum gravity theory to answer how, why and 
what a particular physical constant (like G) takes a specific value. As we are considering only 
variations of G with respect to temporal coordinate t, there is a suggestion to take G(r), as a 
spatial (or radial) dependent quantity [13]. In 14j, it is shown that G was actually varying 
before the electroweak era in the early Universe, while after the spontaneous symmetry 
breaking, the gravitational coupling G attained a finite constant value known today. This 



hue of reasoning can be extended to other types of couplings in particle physics as well M 
but we shall not go into that. In the same manner, some authors [16j have also considered 
the speed of light as a time varying quantity. 

There have been many proposals in the literature attempting to theoretically justified 
a varying gravitational constant, despite the lack of a full, underlying quantum gravity 
theory. Starting with the simple but pioneering work of Dirac 17|], the varying behavior in 
Kaluza-Klein theory was associated with a scalar field appearing in the metric component 



corresponding to the 5-th dimension [l8 and its size variation 



19| . An alternative approach 



arises from Brans-Dicke framework 20(], where the gravitational constant is replaced by 
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a scalar field coupling to gravity through a new parameter, and it has been generalized 



to various forms of scalar-tensor theories 



21], leading to a considerably broader range of 



variable-G theories. In addition, justification of a varying Newton's constant has been 



established with the use of conformal invariance and its induced local transformations 



22]. 



Finally, a varying G can arise perturbatively through a semiclassical treatment of Hilbert- 
Einstein action [23|, non-perturbatively through quantum-gravitational approaches within 
the "Hilbert-Einstein truncation" j^j], or through gravitational holography 25, 26]. 

There is convincing evidence of a varying Newton's constant G: Observations of Hulse- 
Taylor binary pulsar B1913 + 16 gives a following estimate < G/G ~ 2 ± 4 x 10~ 12 yr~ l 



271 ] . helioseismological data gives the bound < G/G ~ 1.6 x lti~ l2 yr~ l [28| (see Ref [29] 
for various bounds on G/G from observational data). The variability in G results in the 
emission of gravitational waves. In another approach, it is shown that G can be oscillatory 



with time 



301 ] . It is recently proposed that variable cosmic constants are coupled to each 



other i.e. variation in one leads to changes in others 31] . A variable gravitational constant 
also explains the dark matter problem as well {32]. Also discrepancies in the value of Hubble 
parameter can be removed with the consideration of variable G 33] . Motivated by the above 
arguments, we shall calculate the the varying G corrections to the statefinder parameters 
for four chosen models of dark energy in the coming section and subsections. 



II. THE MODEL 



34] 



The action of our model is 

S = I d A x£ 



d A x 



-9 



R 

G 



+ F(G) 



(1) 



where G is the Newton's gravitational constant and F(G) is an arbitrary function of G while 
£ m is a matter Lagrangian. From the Euler-Lagrange equation 

<9£ 



— - V 

dG~ "d(0 M G)' 

to obtain 

d£_ B_ 
dG~CP' 

Moreover from the variation with respect to g^, we have from (JTJ) , 



(2) 



(3) 



(4) 
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Writing 



we get from (TJJ 



1 



GF(G) = A(t), 



Sahni et al 



Riiv - ^QtivR = 87rG(t)T M1/ + ^A(t). 



(5) 



(6) 



35| introduced a pair of cosmological diagnostic pair {r, s} which they termed 



as Statefmder. The two parameters are dimensionless and are geometrical since they are 
derived from the cosmic scale factor alone, though one can rewrite them in terms of the 
parameters of dark energy and matter. Additionally, the pair gives information about dark 
energy in a model independent way i.e. it categorizes dark energy in the context of back- 
ground geometry only which is not dependent on the theory of gravity. Hence geometrical 
variables are universal. Also this pair generalizes the well-known geometrical parameters 
like the Hubble parameter and the deceleration parameter. This pair is algebraically related 
to the equation of state of dark energy and its first time derivative. 

The statefmder parameters were introduced to characterize primarily flat universe (k = 0) 
models with cold dark matter (dust) and dark energy. They were defined as 

a 



r = 



a# 3 ' 
r — 1 



3(g 



(7) 
(8) 



Here q = —-^jp is the deceleration parameter. 

For cosmological constant with a fixed equation of state (w = —1) and a fixed New- 
ton's gravitational constant, we have {1,0}. Moreover {1,1} represents the standard cold 
dark matter model containing no radiation while Einstein static universe corresponds to 
{oo, — oo} [3^]. In literature, the diagnostic pair is analyzed for various dark energy candi- 
dates including holographic dark energy 37], agegraphic dark energy [38], quintessence 39], 



dilaton dark energy 40J, Yang-Mills dark energy 4lj, viscous dark energy 42 



interacting 



dark energy 48], tachyon 43], modified Chaplygin gas 44J and f(R) gravity 45( to name a 
few. 

For the present homogeneous, isotropic and spatially flat universe containing dark en- 
ergy and dark matter (ignoring baryonic matter, radiation, neutrinos etc), the Friedmann 
equation is 



H 2 



8nG 



(Pm + Px) 



(9) 
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Let us first consider dark energy obeying; an equation of state of the form p x = wp x . The 
formalism of Sahni and coworkers [35] will be generalized to permit varying gravitational 
constant. In this case, the definition of s is generalized to 

r — Q , . 

ss i^fy (10) 

Here Q = Q m + Q x . The deceleration parameter may be expressed as 

q = ±[n m + (i + 3w)n x }. (ii) 

Hence if Q x , Q and q are determined by measurements the equation of state factor w may 
be found from 



Differentiation of parameter q, together with (JJJ) leads to 



w=-±——. (12) 



From ffTTl). we have 



Furthermore, 



with 



r = 2q 2 + q-j ! . (13) 



1-1 -3 

q = -Q m + -(1 + 3w)Q x + -wQ x . (14) 



^ = — ^-Pcr, (15) 



Per P. 



ci- 



and 



Hence 



.ff G \ 

Per = Per [2— ~ -J, (16) 



ij = -# 2 (l + g). (17) 



p a = -Hp a [2(l + q) + A G ], (18) 
where A G = G'/G, G = HG' which leads to 

Q = ^ + QH[2(1 + q) + A G }. (19) 

Per 

For cold dark matter p m = —3Hp m , (fl9~j) gives 

f2 m = fi m i/(-l + 2g + A G ), (20) 
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and for dark energy p x = — 3(1 + w)Hp x , ( fl9|) gives 



tt x = n x H(-l -3w + 2q + A 



G 



(21) 



Inserting (1201) and (1211) into (1141) and the resulting expression into (1131) finally leads to 



T ~~\~ 



9 

1 + -w(l + w) 



A 



wVt x — [fi m + (1 + 3w)f2 3 



2# 



Inserting the expression ( 122]) into fflUj) leads to 



s = 1 + w — 



w 



A 



G 



[n m + (i + 3w)o 2 



3wi/ 9wVl x 

Assuming a spatially flat universe satisfying Q m = 1 — Q x , we obtain from 
the following expressions 

9 3 Aq 
r = 1 + -tu(l + - -^jjWVL x — (1 + 3wtt x ), 



(22) 



(23) 



and 



s = 1 + w 



w 



A 



G 



3wH 9w^l, 



(24) 



(25) 



Note that and (|25|) provide necessary variable G corrections to statefinder parameters 



of Sahni et al 



351 ] . In particular for the ACDM model, the above pair gives 

A r 

r = — 1-30, , 



90/ 7 



(26) 



(27) 



Moreover the standard cold dark matter corresponds to { — Aq/2, — oo}. For a universe 
containing only dark energy i.e. Q m = 0, we get from (122]) and (123]) : 



9 

1 + + iw) 



3 A 
O z - 2H™Qx 2~Q + 3w )^ 



s = 1 + w — 



A 



G , 



3wH 9w 



l + 3w) 



(28) 
(29) 



A. Interacting dark energy 

There is a class of cosmological models in which evolution of the universe depends on the 
interaction of cosmic components like dark energy and dark matter (and possibly radiation 
and neutrinos). Models in which the main energy components do not evolve separately but 
interact with each other bear a special interest since they may alleviate or even solve the 
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"cosmic coincidence problem" . The problem can be summarily stated as "why now?" , that 
is to say: "Why the energy densities of the two main components happen to be of the same 
order today?" 

The conservation equations for the interacting dark energy-dark matter are 

p x + 3Hp x (l + w) = -Q, 

p m + 3Hp m = Q. (30) 

The coupling between dark components could be a major issue to be confronted in studying 
the physics of dark energy. However, so long as the nature of these two components remain 
unknown it will not be possible to derive the precise form of the interaction from first 
principles. Therefore, one has to assume a specific coupling from the outset 46|, |47J or 
determine it from phenomenological requirements 48]. Here we take Q = — 3UH which 
measures the strength of the interaction j^sj. For later convenience we will write it as 
Q = — 3UH where the new quantity n has the dimension of a pressure. 
The second Friedmann equation is 

H = -A7rG(p t +p), (31) 

where p t = p m + p x . Notice that if the cosmic energy density and pressure violate the null 
energy condition pt+p < 0, then H > 0, i.e. the Hubble horizon if -1 will be shrinking. This 
situation can arise in a phantom- type dark energy model. Also if p t + p > then the Hubble 
horizon is expanding and this can be due to a quintessence-like dark energy. Moreover for 
a cosmological constant dominated universe {pt + P = 0) or during inflation in the early 
universe, the Hubble horizon is a fixed quantity. Differentiating (13T1) w.r.t t and using ( 1301 
and (jSJ), we get 



H 9 / p \ 9 r ,„ . p x tull wp x i '6 . /„ p \ 

— = - l + w (l + w ^ -p- - -A G ! + -)■ 32 

H 4 2\ p t J 2L p t p t 3Hp t \ 2 V p t J 



3 A P 

Pt' 

At variance with H and H, the second derivative H does depend on the interaction between 
components. Consequently, to discriminate between models with different interactions or be- 
tween interacting and non-interacting models, it is desirable to characterize the cosmological 
dynamics additionally by parameters that depend on H. 
Consequently r can be written in alternative form as 

r = j^-3q + 2. (33) 
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Using (1321) in (133)) and the expression of q 



■| — I 2 -), we obtain 



9 w 

1 + ~: 



W7 



n 



2 G V 1 



(34) 



21 + kA p x 3Hw/ 2 \ i + /t, 

where k = p m /Px is a dimensionless but a varying quantity. Note that a constant k require 
extreme fine-tuning. Moreover for a fixed k, the cosmic coincidence problem is trivially 
resolved since the universe will remain in that state forever. To resolve the coincidence 
problem, the density ratio k has to be slowly varying over a time of the order H^ 1 . 
Making use of ( 134)) in (jH)), we obtain 

s = 1+u, -£-d^-tK 1+K+ 4 (35) 

Note that for A G = 0, the above expressions f )34|) and (135)) reduce to the expressions studied 



m 



48| . therefore ( 134")) and (135)) provide necessary variable G corrections to the statefinder 



parameters for an interacting dark energy model. 

In a paper 49j], the authors showed that scaling solutions of the form k ~ a~^, where 
£ denotes a constant parameter in the range [0, 3] can be obtained when the dark energy 
component decays into the pressureless matter. These solutions are interesting because they 
alleviate the coincidence problem 50j. Indeed a model with £ = 3 amounts to the ACDM 
model with w — — 1 and II = 0. As mentioned earlier, the £ = trivially resolves the 
coincidence problem, which is of no interest. Any solution £ < 3 renders the coincidence 
problem less acute. In that scheme, with w =constant, it has been shown in 49] that the 
interactions which produce scaling solutions are given by 



n = / £\ k (i + z)* 

Px r + 3)i + K (i+z)t : 



(36) 



where n = k (1 + z)Z, z = (ao/a) — 1 is the redshift and = p m o/p x o is the present density 
ratio. Putting ( 136)) in ( 134)) and (135)) . we get 



w 



21 + k (1 + z)Z 

3 = 1 + W - 



1 + W 



w + — 



w + - 



K (l 



3/ 1 + K (l + z)tl 

K {1 + zf A G 
3J 1 + k (1 + 3w 



3 A 

-A G 
2 G 



1 + - 



1 + K (l + Z)t 
I + W + K (l + Z)t 



(37) 
(38) 



B. Holographic dark energy 



The holographic dark energy is constructed in the light of the holographic principle. Its 
framework is the black hole thermodynamics 5l| and the connection (known from AdS / CFT 
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correspondence) of the UV cut-of of a quantum field theory, which gives rise to the vacuum 



52j . Thus, determining an appropriate 



energy, with the largest distance of the theory 
quantity L to serve as an IR cut-off, imposing the constraint that the total vacuum energy 
in the corresponding maximum volume must not be greater than the mass of a black hole 
of the same size, and saturating the inequality, one identifies the acquired vacuum energy 
as holographic dark energy: 

3c 2 1 

P* = ^T¥- (39) 



8nG L 2 ' 

Here c is the holographic parameter of order unity. Note that the IR cut-off L can be 
of several types such as particle horizon, Hubble horizon and future event horizon. The 
definition of holographic dark energy is sensitive to the choice of each horizon and gives 
different dynamical features corresponding to each horizon. For the present study, we don't 
need to chose a specific length scale and work with general L. 
The time evolution of (13"9"|) is 

2 



Pol 



(40) 



Using (1401) in the energy conservation equation yields 



1 

w = - 
3 



i - -Jn 



A 



G 



(41) 



Differentiating (14T1) . we obtain 



dw 



H d\\ia 
where we have used 



3c 



— v/fi T - A 



G 



1a; 



Q, r - A 



Putting (S2D and (Ii3"j) in (1211) and (j^, we obtain 
r = 1 + ^-1-^v^+A g )(i-^v^+^g) +^ /2 (l-fi,)(l 



+-A' Q x - — 
3 2 



l + nj - i - - va 

c 



A, 



(42) 
(43) 

-A G 

(44) 



2 
3 

x 



/a 



3c 



+ 



A 



G 



+ 



a; 



G 



1 + 



- 1 | A G 



1 + - 1- 



fi x + A 



G 



a. - a 



G" 



(45) 
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C. New-agegraphic dark energy 



An interesting attempt for probing the nature of dark energy (DE) is the so-called "age- 
graphic DE" (ADE). This model was recently proposed [53] to explain the acceleration of 
the universe expansion within the framework of a fundamental theory such as quantum grav- 
ity. The ADE model assumes that the observed DE comes from the spacetime and matter 
field fluctuations in the universe. Following the line of quantum fluctuations of spacetime, 



Karolyhazy et al. 54] discussed that the distance t in Minkowski spacetime cannot be 



known to a better accuracy than St = fit 2 p /3 t 1/3 where is a dimensionless constant of or- 
der unity. Based on Karolyhazy relation and Maziashvili arguments js5|, Cai proposed the 
original ADE model to explain the acceleration of the universe expansion [53J]. Since the 
original ADE model suffers from the difficulty to describe the matter-dominated epoch, a 
new model of ADE was proposed by Wei and Cai 
be the conformal time n instead of the age of the universe. 
The definition of NADE is given by 56 1 

3n 2 1 



561 ]. while the time scale was chosen to 



Px SvrGV' 

where n is a constant of order unity. Its time evolution is 

2 

Px 



-HpJ—^+Ag). 

\3na J 



Making use of ( 14 7p in the energy conservation equation yields 

2 



to 



-1 + 



3na 



Ac 
3 



(46) 



(47) 



(48) 



Differentiating (1481) w.r.t. t gives 



H W 



3na 



3na 



f2 x 



A 



G 



3na 



tlx + -A 



G- 



where we have used 



Using Q3SD and (SSJ) in (EI| and (ES]), we get 
( 2 



—\Al 

na 



A, 



0- 



(49) 



(50) 



9, 

2' 



1 + 



+— nl /2 

na 



3na 
1 



3 / V 3na 



A 



G 



2na 



^ /2 (l-^)(3 



3na 



Ac 
2 



i + n x (A G -3) + — n 3 J 2 



na 
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A r 1 

jq x + — - - 

3na 3 3 

1 



1 + 



3na 



3na 



V x + -A' G 



A 



A G 

0* + — 
3 

G 



3na 



3na 



3Q X (A G - 3) + £fi 3/2 



1 + 0,(A G - 3) + — n 3 J 2 
na 



(52) 



D. Generalized Chaplygin gas 



Interest in generalized Chaplygin gas (GCG) arose when it appeared that it gives a unified 
picture of dark energy and dark matter. It solely provides the density evolution of matter 



at high redshifts and dark energy at low redshifts 
it explains the recent phantom divide crossing 



58j |. Other successes of GCG is that 



supernova 



59], is consistent with the data of type la 



601 ] and the cosmic microwave background [61]]. The GCG emerges as an effective 



fluid associated with d-branes 



62] and can also be obtained from the Born-Infeld action 63 ]. 



In the present context, we are treating GCG as a dark energy candidate. 
The equation of state for generalized Chaplygin gas is 

A 

Px — ~J 
Px 

where a and A are constants. Using ( 1531) in the energy conservation equation yields 



(53) 



A + 



B \ l+c 



a 3(l+a) 

Here B is an integration constant. The corresponding state parameter takes the form 

B 



(54) 



w 



-AIA + 



a 3(l+a) 



(55) 



whose time derivative yields 



w 



-3ABH(l + a)[A + 



B 



a 



3(l+a) 



(56) 



Using 055]) and ([56} in fl24l and rt25J, we get 



r = l + 9 -ABn x (l-a-^)(A+-^ a ;- 2 



Ag 
2 



1 - 3AVL X A + 



B 



a 3(l+a) 



.(57) 



i -^ A+B ^ + < A+ j^y + ^-X A 



B 



a; 



3(l+a) 



1 -3Q x a(a + 



B 



a 



3(l+a) 



-In 



(5f 
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III. CONCLUDING REMARKS 

In this paper, we calculated the corrections to statefinder parameters due to variable 
gravitational constant. These corrections are relevant because several astronomical obser- 
vations provide constraints on the variability of G. An important thing to note is that the 
G— corrected statefinder parameters are still geometrical since the parameter A G is a pure 
number and is independent of the geometry. 
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